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SUMMARY • 


A theoretical investigation h£is been made to determine the effects 
of an elastic -wing on the dynam i c longitudinal stability of thin-iri.ng 
airplane configurations. In order to investigate the effects of various 
important parameters, the configurations were assumed to vary in wing 
sweep angle from 0 ° to 6o°, in center-of- gravity location from 25 per- 
cent mean aerodynamic chord to L5 percent mean aerodynamic chord, and 
in ratio of iring mass to airplane mass from O.I5 to O.5O. 

Three degrees of freedom were assumed - freedom in vertical trans- 
lation of the rigid airplane, pitching rotation of the rigid airplane, 
and displacement of the wing tip due to bending of the elastic wing. 

The elastic wing mode was determined from consideration both of the 
deflection under static loading and the deflection in the primary ground- 
vibration mode and was represented by a combination of bending in the 
primary mode and the associated torsion. Lagrange's method was -used 
to obtain an equation of motion for each degree of freedom. The char- 
acteristic equation of the system was solved for the period and damping 
of the airplane mode and of the wing mode. Solutions were also obtained 
for three sin^jllfied airplane and wing systems: (l) the wing mode alone, 

(2) the airplane under quasi-static conditions, and (5) the rigid airplane. 


An analysis of the solutions showed that, for configurations having 
to 60 sweepback, no dynamic instability due to wing flexibility was 
indicated; however, the loss in static stability due to wing flexibility 
was found to be a fairly serious problem. For configurations having no 
sweepback, the wing was subject to a decrease in oscillatory stability 
for the large ratio of wing mass to airplane mass accon5)anled by forward 
center-of-gravity locations. The quasi-static method gave results cont- 
peirable to those of the semirigid method for sweptback wings; however, 
for straight brings the quasi-static method gave poor results. 
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INTROrOGTION 


As a resiilt of the trend towaird the use of thin swepthack wings on 
large hi^-speed aircraft, the effect of structural flexihility on dyna- 
mic longitudinal stahiUty must he considered. The inclusion of flexi- 
hility modes as degrees of freedom in addition to those of the rigid 
airplane results in hi^er order characteristic equations. As the order 
of the characteristic equation increases, the amount and difficulty of 
work involved in obtaining the solution increases . The additional work ■ 
is, of course, unjustified if the results do not differ appreciably from 
the results obtained by nsin g less rigorous methods. 

The present paper presents a theoretical investigation of the 
effects of wing flexibility on the dynamic longitudinal stability of 
airplanes by including a degree of freedom for the elastic wing. The 
resTolts of the higher order equations are con5)ared with the results of 
less rigoroois methods. In oitLer to investigate the effects of veirious 
in^jortant paraneters, the configurations are assumed to vary in wing 
sweep angle from 0 ° to 6 o°, in center- of- gravity locations from 25 per- 
cent mean aerodynamic chord to 45 percent mean aerodynamic chord, and 
in ratio of wing mass to airplane mass from O.15 to O.5O. A single 
value of wing bending of torsional stiffness typical of those pos- 
sessed by an actual airplane is used in the analysis . 

The equations of motion are derived in the appendix by using 
Lagrange's method. The resulting semirigid three-degree-of- freedom 
equations of motion are solved at various fli^t conditions for the 
period and damping of the two modes of oscillation, the first involving 
primarily wing motion and the second, airplane motion. These solutions 
are con^jared with solutions obtained for three slii5)lified airplane and 
wing systems: (l) the wing mode alone, (2) the airplane tinder quasi- 

static conditions (flexible inertia and danplng terms assumed to be 
zero, and (5) the rigid airplane. 


SYMBOLS 


A <11 pgrnTTi showing the system of axes and positive directions of 
forces anti moments on the airplane is presented in figure 1 . 

(b cos a)^ 

A aspect ratio, 


generalized nondlmensional mass coupling term between Z and 
h degrees of freedom, a^t^/ pSc 
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generalized nondlmens lonal mass term of flexible wing mode 
between elastic wing and b degree of freedom, a^pSc 

generalized nondimensional mass coupling term between 0 and 
b degrees of freedom, a^-^/pS^ 

generalized mass coupling term between Z and b degrees of 


To 

freedom, 2 J'^ [m^'f^Cy) - S^'f^Cy^dy^, slugs 


generalized mass term of flexible wing mode between elastic 
wing and b degree of freedom, 

slugs 

generalized mass coupling term between 0 and b degrees of 
freedom, 

lw*^0(y) - “w'ef2(y) + S^'ef0(y^dyo, 

slug- ft 

span (along elastic axis), ft 

force coefficient due to elastic wing deflection, Fj^/qS 
nor m al-force coefficient, N/qS 
trim normal- force coefficient, W/qS 

pitcblng-moment coefficient about Y-axis, M/qSc 
mean aerodynamic chord, ft 
chord, ft 

section pitching-moment coefficient about wing elastic axis, 

m'/q.c^ 

section normal-force coefficient, n/qc ■ 
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■bending stiffness, lb-in. ^ 
kinetic energy, ft-lb 
potential energy, ft-lb 

longitudinal distance from airplane center of gravity to wing 
elastic axis (function of spanwlse location), positive 
forward, ft 

force , lb 

section force, lb /ft 

spanwlse bending mode shape along wing elastic axis 

spanwlse twisting mode shape about wing elastic axis per unit 
tip bending deflection, radians /ft 

torsional stiffness, lb-in. ^ 

acceleration due to gravity, ft/sec^ 

wing-tip deflection, h/c, chords 

wing-tip deflection of elastic axis due to bending, positive 
downward, ft 

moment of inertia about Y-axis, m^x 2 , m^xS, or nv^ 2 ^ 
slug-ft^ 

section moment of inertia, m^'x^^ slug-ft^/ft 
radius of gyration about Y-axis, chords 
reduced angular frequency, (oc/'Y 

longitudinal distance from quarter chord of wing mean aero- 
dynamic chord to quarter chord of tail mean aerodynamic 
chord, ft 

pitching moment about Y-axis, ft-lb 

section pitching moment about Y-axis, ft-lb/ft 

mass, slugs 

section pitching moment about elastic axis, ft-lb/ft 
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section mass, slugs/ft 

N normal force, positive downward, Ib 

n section normal force, Ib/ft 

Q generalized coordinate 

g dynamic pressure, Ib/sq. ftj also, pitching angular velocity, 

radians/sec 

rg ratio of local chord to root chord, ■ c/cr 


S 

s' 


d 


^A^Sf ,Sy^ 


t 

V 

W 

X 


X 

Y 

7 

Z 


z 


a 

6 


wing plan-form area, sq ft 

section mass moment about elastic axis, ™w'x, slug-ft/ft 

mass moment, m^x, or i%x, slug-ft/ft 

time to danp to 0.1 amplitude, sec 

time, sec 

velocity, fps 

weight, lb 

longitudinal axis of displacement fixed at airplane center 
of gravity 

longitudinal displacement, positive forward, ft 

lateral axis of reference fixed at airplane center of gravity 

/ 

lateral or spanwlse displacement, ft 

vertical displacement of airplane center of gravity, positive 
downward, ft 

vertical wing deflection of elastic axis due to wing 
bending, positive downward, ft 

angle of attack, positive wing leading edge up, radians 

elevator deflection, positive trailing edge down, radians 

dimensionless spanwise coordinate, fraction of semispan 
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dimensionless spanwise coordinate, — fraction of exposed 


semlspan. 


^0/2^ 


A 

V- 

P 

cu 


angle of pitch shout airplane center of gravity, positive 
airplane nose up, radians 

sweep angle of eleistic axis of wing, deg 

nondimensional airplane mass, m^/pSc 

mass density of air, slugs/cu ft 

angle of twist of airfoil in plane perpendicular to elastic 
axis, positive wing leading edge up, radians 


angular frequency, radians /sec 


Subscripts: 

A airplane 

av average 

eg center of gravity 

f fuselage 

h flexihle-wlng degree of freedom 

i intersection of elastic axis with fuselage 

max TnavT miTm 

o exposed wing 

Q generalized coordinate 

r wing root 

t tail 

w wing 

Z vertical degree of freedom 

0 pitching degree of freedom 
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Dots are used to indicate differentiation with respect to time; 
• d0 

for example, 9 = — . 

dt 


The subscripts a, 0, h, h, q, H, and 6 indicate differentia- 

dCjT 

tion with respect to the subscript; for example, Ctj = . 

a, da 


BASIS OF ANAIZSIS 
Semirigid Msthod 


A flexible structure may he considered to have an infinite number 
of degrees of freedom because it can deflect an infinite number of ways 
depending on the loading. However, in order to make the solution for 
the motion of the flexible airplane as simple as possible, the number 
of degrees of freedom and consequently the nunDber of equations of motion 
should he kept to a minimum. This method of approach is presented in 
reference 1 and is called the semirigid concept. In this approach the 
mode of flexure under load is always the same regardless of the loading. 

Three degrees of freedom were assumed for the configuration studied 
herein - freedom in vertical translation of the rigid airplane Z, 
pitching rotation of the rigid airplane 0, and displacement of the wing 
tip due to bending of the elastic wing h - and an equation of motion 
must he derived for each of these three degrees of freedom. A slnplified 
way of deriving these equations is hy Lagrange's method which is described 
in reference 2. Lagrange's method consists in writing an ejqoression for 
the total kinetic energy and the total potential energy of the system. 

The operations indicated hy the general Lagrangian equation are then 
performed in the expression for total energy, and an equation of motion 
for each degree of freedom is obtained. The application of this method 
to a configuration with a flexible wing is presented in reference 3* 

For the convenience of those who- are not familiar with Lagrange ' s method, 
a derivation of the equations of motion used in this paper is presented 
in the appendix. 

The flexible wing mode shape was assumed to consist of bending fz(y)^ 
combined with twisting per unit bending deflection at the wing tip f^^y) . 
When the wing is flexed in the degree of freedom characterized by the^ 
wing-tip bending deflection h, the bending deflection at station y is 
z and the torsional deflection is 0 . Ho motion of the fuseleige was 
considered to be associated with the wing deflection mode. In prantice, 
some fuselage rotation or bending woxild occur in conjunction with the 
wing mode and should be considered when making a specific analysis. 
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In order to use the equations of motion, the characteristics of a 
wing mode which satisfactorily represents the flexing configuration must 
he determined, and the conditions under which the wing flexes must he 
specified. The wing was assumed to he a flexible cantilever beam which 
was anchored to the fuselage in a plane through the point of intersection 
of the elastic axis and the fuselage and perpendicular to the elastic 
axis. The basic mode of oscillation was considered to be a conhination 
of primary bending in a vertical plane containing the elastic axis and 
the associated twisting about the elastic axis. The shape of the bending 
and twisting modes couM be expected to vary somewhat with the span load 
distribution on the wing. At zero frequency the loading is static, the 
result of steady additional lift forces. As the frequency increases 
toward the wing natiiraJ. frequency, the wing mode shape would be expected 
to approach the mode shape corresponding to an elastic vibration of the 
wing on the ground. Simple beam theory was used to obtain the static 
bending mode shape and the twisting mode shape under steady load. The 
moment of the additional lift forces about the wing elastic axis was 
used in determining the twisting mode shape and the amount of twist per 
unit tip deflection. The twisting mode was determined solely for static 
loading inasmuch as the twist occtirring in an elastic vibration of the 
wing appeared to be negligible. Th is obsei^ration resulted from the fact 
that, for the assumed wing mass distribution, the wing was almost per- 
fectly mass balanced about the elastic axis. 


The bending mode shape of the wing oscillating at the ground natural 
frequency was obtained with the aid of reference The two calculated 
bending modes and the calculated twisting mode are presented in figure 2. 
The calciilated bending modes closely approximate each other despite the 
large variation in distribution between the static and dynamic loading. 

A parabolic variation of spanwise bending- fz(y) = and a linear 


variation of spanwise twist 


^ no 


are also presented in fig- 


ure 2. The agreement between the calcxilated mode shapes and these sim- 
plified approximating curves is seen to be veiy good; therefore, the 
sin^llfied c\irves were used to represent the mode shapes. 


The values of wing bending and torsional stiffness and the wing 
structural weight distribution for the assumed airplane configurations 
were based on those possessed by an actxial swept-wing bomber airplane 
and are presented in figure 3- The bending and torsional stiffnesses 

vary along the span in accordance with r^^, the fourth poorer of the 
ratio of the local chord to the root chord, as suggested in reference 5» 
This assumption is in good agreement "i-Jlth the structural characteristics 
of the actual airplane. The actual aiiplane had a ratio of wing struc- 
txiral mass to airplane mass of O.I 5 . The higher values of the ratio of 
wing mass to airplane mass of 0.55 and O. 5 O as presented in figure 5 
approximate some typical wing-airplane mass ratios for wings having 
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ELddl 131011811 maBses attached to them such as nacelles and exrtemal or 
internal stores. However, in these cases it should he noted that the 
eissumed mass distribution does not necessarily correspond to an airplane 
having such additional masses attached to it. 

The bending mode shape was used to calculate the ground natural 
frequency of the wing. In order to obtain a range of natural frequencies 
of the wing, the wing stiffness was assumed to remain constant nnil the 
mass of the wing was assumed to increase so that the airplane mass also 
increased. Frequencies were calculated for ratios of wing mass to air- 
plane mass of 0.15, 0 . 33 } and O.5O and are presented in table I. 

Additional configurations used in the analysis were obtained by 
modifying a beisic configuration. The variations were obtained by holding 
the wing area and the span along the elastic axis constant as the wing 
panels were swept about a vertical axis which was located at the Juncture 
of the wing elastic axis and the airplane center line. As the wing was 
varied in sweep, the wing root was moved forward or rearward so as to 
keep the tall length l-j; constant. Configurations having sweep angles 
of 0°, 40°, and 6o° measvired with respect to the elastic axis were assumed 
and are shown in figure 4. Because of the large aspect ratio, only a 
small difference in angle of sweep between the elastic axis and the 
quarter-chord line exists. In view of this small difference anti for pur- 
poses of convenience, the angles of sweepback used herein refer to the 
sweep of the elastic axis which is the 38-percent-chord line. The fuse- 
lage and tail dimensions were held constant. The airplane center of 
gravity was assumed to be located at 25, 35, and 4-5 percent of the wing 
mean aerodynamic chord. 

The inertia terms and the aerodynamic dangping and restoring-spring 
parameters were calculated for the configurations and are presented in 
table II. The airplane lift and pitching-moment characteristics were 
calculated by using references 6 and 7. Strip theory tmcorrected for 
compressibility effects was used in calculating the force parameters . 

In the con5)utatlons of the generalized mass terms, some small terms were 
neglected. The equations of motion were solved for period nnil dan5)ing, 
and solutions for each configuration were obtained at three or more 
dynamic press\ires. Also, transient solutions for a configuration at 
two altitudes were presented in order to show the relation of the wing 
oscillation to the airplane oscillation. Because of the variations with 
altitude of the relative density of the airplane and the generalized 
meisses , the solutions of the characteristic equation for values- of dyna- 
mic press\n:e at one altitiide are not applicable for the same values of 
dynamic pressure at other altitudes. 
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Other Methods 

In order to determine whether more simple methods could he used to 
predict the dynamic characteristics of the flexible airplane and wing, 
results were obtained by the quasi-static method and the wing-mode -alone 
method. The quasi-static method consisted in eliminating all the inertia 
and damping terms pertaining to the flexible mode (temus containing 
and DH) from the equations of motion and solving the resulting second- 
order characteristic equation for period and damping. In effect, the 
results of this method are comparable to those of the rigid-airplane 
method except that the effects of flexibility on the aerodynamic param- 
eters are accounted for at any given dynamic pressure. It should be 
noted that the rigid-a±rplane solutions are the same as the solutions 
of the semirigid airplane at zero dynamic pressure. The simple wing- 
mode-alone method consists in neglecting the first two equations of 
motion and solving the third equation of motion after neglecting the 
coupling terms (terms containing Da, a, D^, and D0) that appear in 
the third equation. It should be noted at this point that the condi- 
tions for the wing-mode-alone method do not simulate the conditions used 
in flutter work. The consideration of a separate torsional mode (as in 
flutter work) is beyond the scope of the problem dealt with in this paper. 
The frequencies of oscillation of the modes are low enough so that con- 
sideration of unsteady-3J.ft effects is unwarranted. The wing-mode-alone 
method was evolved in order bo imderstand better the effects of short- 
period coupling in the wing mode. This method, therefore, represents 
the wing motion that would occxxr if the airplane center of gravity was 
constrained to move in a strai^t line at all times without pitching. 


RiSUHES AUD DISCUSSION 


The equations of motion were solved for the daagjing and frequency 
of the flexible wing mode and of the airplane mode for various dynamic 
pressures at a standard altitude of 8,000 feet. The solutions for ratios 
of wing mass to airplane mass of O.I 5 , 0.55; and O. 5 O and for center-of- 
gravity locations of 25 , 35, snd percent of the wing mean aerodynamic 
chord at zero sweep angle are presented in figure 5* sweep angles 

of 40° and 6o° the solutions did not change appreciably with increasing 
mass ratio. Therefore the res\ilts for the mass ratio of 0.33 only and 
with center-of-gravity locations of 35# and ^5 percent of the mean 
aerodynamic chord are presented in figure 6. The wing characteristics 
covered in parts (a) of figures 5 and 6 are expressed in terms of actual 
timej whereas the airplane characteristics covered in parts (b) of these 
figures are expressed in terms of nondlmenslonal time. The results are 
presented in different forms inasmuch as each form is believed to enable 
the best interpretation of the data. The wing groimd natural frequency 
is shown in the plots of wing frequency against dynamic pressure in 
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figures 5 (a) stud 6(a). The wing-mode- alone solutions and the solutions 
for the quasi-static conditions of the airplane mode are also shown in 
these figures. 

The results for the airplane mode are presented as a function of 
the period and of the reciprocal of the time to damp to 0.1 an^litude 
where time is expressed nondimensionally as units of distance traveled 
measured in fuselage lengths. Fuselage length, rather than the usual 
unit of wing chord, was used because the fuselage length remained con- 
stant for n.1 1 configurations whereas the wing mean aerodynamic chord 
vsried with the sweepback angle. Solutions for the rigid airplane con- 
figuration presented in this form, are independent of forward speed. At 
zero dynamic pressvure flexibility effects must vanish, and the solution 
for the flexible airplane mode is identical to that of the rigid con- 
figuration. Therefore, the deviation of the solutions from the values 
for the configurations at zero dynamic pressure represents the effect 
of flexibility at any particular dynamic pressure at an altitude of 
8,000 feet. 

Before examining the data, it is well to note a relation between 
the wing mode and the airplane mode. By expanding the nondimens ional 
characteristic equation, the sim of the danping of the airplane mode 
and the wing mode for a particular configuration can be shown to remain 
constant with variation in dynamic pressure. Any changes in wing 
danping, therefore, result in opposite changes in airplane damping. 


Semirigid Method 

The solutions for the configurations having zero sweep obtained by 
the semirigid method are discussed first. The effects of variation of 
wing mass and center- of- gravity location are shown for zero sweep in 
figure 5 * 111 "the lower half of figure 5 (a) ^ "the danping curves do not 

extend to zero dynamic pressure. In the region of zero dynamic pressure, 
the danping curves would approach infinity because of the assumed absence 
of structural danping. 

The effective danping of the wing oscillation decreases as the mass 
of the wing increases. This effect is shown in the danping curves by 
an increase in time to damp to 0.1 anplltude and also in the frequency 
curves by a tendency for the frequency to become constant at higher 
dynamlq pressiures. The danping for the wings of greater mass also shows 
an appreciable variation with center- of- gravity location at the highest 
dynamic pressures. For the center of gravity at the most forward loca- 
tion ( 0 . 25 c), the oscillation tends to become less stable with increase 
in dynamic pressure. This effect is indicated by an increase in time to 
danp to 0.1 anplltude while the frequency becomes approximately constant. 
Solutions are presented at dynamic pressures which would be above the 
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critical Mach mimher of a configuration in order to gain a general idea 
of the stahillty trends. Although the solutions at those dynamic pres- 
stires are meaningless for the present configuration, they would apply 
at lower dynamic pressures for a configuration with a wing of greater 
flexibility. 

Some effects of the airplane mode of oscillation on the wing mode 
may he seen from a con^axlson of the wing-alone solutions with the semi- 
rigid solutions for the wing. To analyze these effects better the fre- 
quencies of the airplane mode have been put in dimensional form and are 
plotted in figure 5(a) together with the wing frequencies. The airplane 
motion generally decreases the dan^ping of the wing oscillation especially 
for the = 0.50 configuration at the higher dynamic pressures. 

The frequency of the airplane oscillation approaches the nat\aral. fre- 
quency of the wing oscillation at these dynamic pressures. The ten^jo- 
rary increase in wing frequency above the ground natural frequency of 
the wing at low dynamic pressures also may he attributed to the effect 
of the airplane oscillation on the wing oscillation. This conclusion 
is substantiated by the lack of increase in frequency for the wing-mode- 
alone method. As is well-known, danping in a slngle-degree-of- freedom 
system, whether positive or negative, will decrease the frequency from 
that of the system with no damping. 

For the airplane mode (fig. 5(^))^ danplng shows a slight gen- 
eral decrease with increase in wing mass and, hence, airplane mass. For 
the light configuration Dv/m^ = O.I 5 with the center of gravity located 
at 45 percent mean aerodynamic chord, the damping is sufficient to cause 
the airplane oscillation to become critically over damped. The danping 
of the airplane oscillation increases slightly with increase in dynamic 
pressure. Generally, tbe danping of the airplane mode for configura- 
tions having zero sweep angle is good. 

At zero sweep angle the effects on the local incidence angle of 
wing flexibility are due to the wing twisting - the bending deflections 
have no effect. As the sweep angle of a wing increases, the effect of 
twist on the local incidence angle gradually decreases and the bending 
conponent of the flexible mode begins to exert a powerful effect. At 
5° of sweep for tbe assumed basic wing, the effect of twist on the local 
wing incidence is neutralized by the effect of bending; therefore, the 
wing is aeroisoclinlc . Above 5° of sweep the bending effects are 
predominant . 

» 

An examination of figure 6(a) shows that the wing frequency Increases 
with increasing dynamic pressure for wings of 40° and 6o° sweepback. The 
trend of the Increase in wing frequency above the wing ground natiural fre- 
quency at low dynamic pressures is the same as for the zero sweep angles; 
for the 40° and 60 ° swept wings, however, the wing frequency continues to 



NACA TN 5251 


13 


increase apparently as a result of increased aerodynamic restoring 
moments wMcli add to the elastic restoring moments of the wing. In gen- 
eral, the data show that the wing mode is always satisfactorily damped 
for the expected flight range of dynamic pressures. The damping appears 
to become poorer at low dynamic pressures, hut these results are pessi- 
mistic because of the ass^uIption of zero structural damping. It is also 
apparent that the center-of-gravlty position has no appreciable effect 
on the characteristics of the wing mode. The damping is slightly less 
for the 60° sweptback configuration than for the k0° sweptback configura- 
tion because the lift-curve slope. is less for the 6o° configuration. 

The characteristics of the airplane mode for the swept-wing con- 
figurations (fig. 6(b)) indicate that the static- stability considerations 
completely overshadow the dynamic stability characteristics in importance. 
At the lowest dynamic pressures for the most rearward center- of- gravity 
location the semirigid configurations possess an oscillation similar to 
that found in the rigid airplanes. With increased dynamic pressures, 
this oscillation changes into two convergences as the maneuvering neutral 
point approaches the center-of-gravity position being considered. With 
still further increase in dynamic pressure, the maneuver point becomes 
coincident with the center of gravity being considered, and one of the 
convergences changes to a divergence. Any further Increase in dynamic 
pressure causes the divergence to be more severe as the maneuvering 
neutral point moves feirther ahead of the center of gravity. For more 
forward locations of the center of gravUy the same sequence occurs 
except at higher dynamic pressures. This analysis is confirmed by the 
data of figure J, which shows the location of the maneuvering neutral 
point with dynamic pressxire obtained by the quasi-static method for 0°, 
40 °, and 60° sweepback angles. 

The top part of figure 6(b) shows that when an oscillation does 
exist the period increeises with increasing dynamic pressTire as would be 
expected from the decrease in static stability with increase in dynamic 
pressure. A comparison between the variation of wing frequency and 
airplane frequency with dynamic pressure shows that these variations 
are in opposite directions so that interaction between these two modes 
appears to be impossible. The loss in static stability with increasing 
dynamic pressure is caused by a forward shift in wing aerodynamic center. 
There is, of course, a compensating factor in the form of decreased 
wing- lift- ciirve slope due to bending, but this factor is less important 
than the forward shift in wing aerodynamiic center. 


Evaluation of Simplified Methods 

Results obtained from the rigid- airplane solutions have already 
been covered in the previous discussion Inasmuch as the results are the 
same as for the flexible eiirplane at zero dynamic pressure. In order 
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to determine whether more single methods could he used to predict the 
dynamic characteristics of the flexihle airplane and wing, results were 
obtained hy the qjoasi- static method and the wing-mode- alone method, and 
these results are shown in figures 5 arwi 6. 

An inspection of figure 5(^) shows that, for the unswept airplane 
configuration at the highest dynamic pressiores, the period and dandling 
can he determined with only a fair degree of accuracy hy the quasi-static 
method. It is noteworthy that the resvilts of the quasi-static method do 
not differ appreciably from the results obtained with the rigid- airplane 
equations of motion. For swept-wing configurations (fig. 6(b)), however, 
the period and daa5)lng are seen to be in somewhat better aigreement with 
those determined by the semirigid method. The major effects of static 
stability are shown by the quasi-static method. The variations in 
damping with dynamic pressure are almost exactly -the same as those 
given by the semirigid method. It therefore appears that the quasi- 
static method can be used to give a good first approximation of the 
dynamic characteristics of a swept-wing airplane having a flexible wing. 

Application of the wing-mode- alone method to the prediction of the 
period and dan5)ing of the flexible wing mode (figs. 5(a) an<i 6(a)) shows 
that for unswept wings the period and damping can be predicted with a 
fair degree of accuracy except in the case of large wing mass ratios with 
far -forward center-of-gravlty positions. For sweptback wings within 
the ranges of parameters covered, the wing-mode- alone method seems to 
give good results. 


Transient Solution 

The period and d amp ing for two modes of oscillation of a series of 
flexible-wing aircraft have been presented and discussed. In order to 
give a better idea of the relative in^iortance of the wing oscillation 
in the total airplane motion, solutions have been obtained on the 
Reeves Electronic Analog Computer for the transient motions following 
a disturbance. The longitudinal transient responses of a 55° swept- 
wing eiircraft configuration to an elevator step input at altitudes of 
8,000 feet and 30,000 feet are presented in figure 8. The configuration 
closely resembles the 40° swept aircraft for flight configuration of 
rn^/m^ = 0.35 in mass, structural, and aerodynamic parameters. The plot 
of wing- tip deflection h shows only a very small excitation of the wing 
mode ■vdiich quickly disappears. It is possible that the wing oscillation 
v mild be excited to a greater degree by a gust or some other form of 
disturbance; however, these records indicate that the wing oscillation 
is not easily excited by xise of the elevator. The airplane oscillation 
is heavily dan®ed at both altitudes. 
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COKCEUSIOWS 


On the basis of an analysis of the effects of wing flexibility on 
configurations varying in sweep angle, center-of-gravlty location, and 
ratio of wing mass to airplane mass, the following conclusions are 
indicated; 

1. For configurations having to 6o° sweephack, no dynamic 
instability due to wing flexibility was indicated; however, the loss 
in static stability due to wing flexibility was found to be a fairly 
serious problem. 

2 . For configurations having no sweephack, the effects of wing 
flexibility were found to be serious only in the case of configurations 
having a large ratio of wing mass to airplane mass and forward center- 
of-gravity locations; for these cases, the analysis indicated that the 
wing was subject to a decrease in oscillatory stability. 

3. The quasi-static method appears to be fairly realistic as a 
means of obtaining a good approximation to the dynamic characteristics 
of a swept-wing configuration. However, the method does not appear to 
be accurate for unswept configurations at high dynamic pressures. 

!<■. A sin^jllfied wing-mode- alone method was found to give good 
results in predicting the characteristics of the flexible wing mode 
of motion for swept-wing configurations. However, this method did not 
appear to be reliable for use with unswept configurations. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 26, 195^ • 
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APPENDIX 

DERIVATION OF ECSJATIONS OF MOTION 
The Lagrangian equation is 

^aQ / 3 q 3q 


(Al) 


where Ejj. and Ep are the kinetic and potential energies of the dynamic 
system. The kinetic energy of the system is equal to the sum of the con- 
tributions of the fuselage and of the elements of the wing which are 
external to the wing-fuselage intersection- In this analysis the fuse- 
lage is assiimed to be rigid, and no motion of the fuselage is considered 
to be associated with the wing deflection mode. The motion of the fuse- 
lage is therefore due to motion of the rigid airplane, and the motion of 
elements of the wing is due to motion of the rigid airplane and motion 
of the flexible wing. The tail aissembly is considered part of the fiise- 
lage. The velocities and the geometric relations of the fuselage and 
wing elements to the airplane center of gravity are presented in fig- 
ure 9- From figure 9(a) the velocity at the center of gravity of the 
fuselage if is seen to be 


if — z — xB (ap) 

and the kinetic energy of the fuselage unit is 

Bi^ = I mf(Z - x8)^ + i I^g^0'2 (A5) 

By expanding equation (A5)^ the kinetic energy of the fuselage becomes 

Ekf = I nifZ2 _ ^ 1 


where 


If0^ = mf(xe)^ + Icg^e 
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For the wing, the vielocity at the center 
element of wing may he obtained with the 


of gravity of a small 
aid of figure 9(^) as 


Zw = Z - (e + x)e + fz(y)h - f^(y)ix (A 5 ) 

The kinetic energy of a RTrm 1 1 element of iring is 

•^k^ = I %'[z - (e + x)0 + fz(y)h - f 0 (y)h^ <3y + | ^cg^^' 2 <3y 


(A$) 


In the present analysis the term 


2 ^cgw 


0 + f0(y^ <3y 


has been con- 


sidered to be small for the wings of high aspect ratio which have been 
considered herein and, therefore, has not been carried further in the 
derivation. Equation (a 6) is integrated over the exposed semispan bo/2. 


To put the equation in proper form for integration, let the variable 
over the exposed semispan be yo-= y - 71 * Tke equation is then inte- 
grated from 0 at the wing- fuselage intersection to bo/2 at the wing 
tip. Expanding equation (a 6) resvilts in 





- 2Ze9 + 2Zfz(y)h + (e0) 


2e0fz(y)h + 


[fzCy)^ + Zf0(y)h - e0^ - e0f0(y)h + 0fz(y)h + 

fz(y)f0(y)tL^ + I Iw'^e^ + 2f0(y)0h + |f0(y)^ |)dyQ (AT) 


Equation (A 7 ) contains mass terms pertaining to the rigid and 
flexible degrees of freedom. These relative mass terms can be grouped 
and assigned definitions. The terms of equation (A 7 ) pertaining to the 
rigid degrees of freedom are coi^ined with equation (A^J-) to obtain the 
kinetic-energy equation of the rigid airplane. The mass terms in equa- 
tion (A 7 ) pertaining to the rigid degrees of freedom are obtained by 
removing the mass terms containing the flexible degree of freedom £ 
from the equation. Then the kinetic energy of the rigid wing is obtained 
an 
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= 2 J ° . I iVi^ - 2e0Z + (e0)5- Sw'(Z0 - e0^) + | 

(a8) 

If the terms of equations (a?!-) and (a 8) are combined, the kinetic 

energy of the rigid airplane in vertical translation hecomes 


I = I “fZ^ ^ ^ fo° i 


(A9) 


Similarly, if the 0^ terms of eqaatlons (Ah) and (a 8) are combined, 
the kinetic energy of the rigid airplane in pitch hecomes 


1 T _ 1 t -^2 

5 Ia® - 2 ^ ^ 


f I m„’(e§)^ + Sw’e02 + ^ dy^ (AlO) 


The terms remaining in eqaatlons (Ah) and (a 8) can he recognized as 
static unbalance terms. Because eqaatlons (ih) and (a 8) are derived 
about the airplane center of gravity, the sum of the static unbalance 
terms must he equal to zero or 

-SfZ0 - 2 f (nv'eOZ + S,^’:M)(ayo = 0 


Combining eqaatlons (A9) and (AIO) yields the total kinetic energy of 
the rigid airplane in vertical translation and pitch 


Ek = I 



(All) 


Tn a similar manner, some nev terms representing the Inertia effects 
of the fleyihle airplane may he defined. After the terms in equations (AT) 
containing ir are collected, the generalized mass pertaining to the 
flexible degree of freedom may he defined as 


•hh^^ = 2h^^ ° [fz(y)I|^ - 2Sw'fz(y)f0(y) + iv' [f0(y^ ^|'3yo 


(A12) 
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By coll acting the mass terms from equation (A 7 ) containing ^ and 0h, 
the generalized masses pertaining to the coupling terms may likewise he 
defined as 


2a7hZh = 2Zh< 


tvr'^z(y) - Svr'f0(y^dyo| 


(AI 3 ) 


2aQj^0h = 20*h ja J' [^S^’fgCy) + I,f'f0(y^dyo - 




(Al4) 


The ahhreviated equation for Mnetic energy is then 

= iSLp^?' + Ti.^^ + a>i>ih + 2a^-^Zh + 2a0j^0h (AI 5 ) 

The potential energy of the system is congiosed of a contribution 
from the airplane and a contribution from the flexible wing. The 
potential energy of the aiip)lane due to its vertical position is given 
by -ZW^. The potential energy of the flexible wing is expressed in 
terms of the frequency of the wing oscillating in the assumed deflec- ' 
tion mode. If the wing is performing a sinusoidal oscillation, then 
the potential energy at. its point of maximum deflection is equal to 
the kinetic energy of the wing as it passes through the point of zero 
deflection. For a small strip of wing vibrating sinusoidal 1y in the 
deflection mode, the maximum velocity is coCzmnv - 0max^^ 
corresponding potential energy expressed in terms of kinetic energy is 

2 

" 2 [^(^ijiax " ^max^^ ^o 


For any spanwise location, z^ny = f^Cy)! and = f0(y)h; substi- 

tuting these values and integrating yields 




1 

2 



“v' [£z(y[l^ - 2Sw’fz(y)f0(y) 



= I oAiSaijii 
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If the potential energy of the airplane is conibined with the potential 
energy of the elastic wing, the total potential energy of the system 
can he expressed as 

Ep = -m^gZ + ^ (AI 6 ) 

where = m^g. 


The equations of motion for each degree of freedom determined hy 
substituting equations (AI 5 ) and (A16) into the Lagrangian equation (a 1) 
and using the genera l ized coordinates Z, 0, and h are 


®Zh^ “ " ^Z 

(ait) 

S0h^ = ^0 

(AI 8 ) 

s-hhii + azh^* + ^Oh® + ahh^li = ^h 

(AI 9 ) 


The generalized forces ^ 0 ) ^h a.ccount for the forces not 

included in the potential-energy term Ep. In order to obtain the com- 
ponents of the generalized force term, the airplane is assumed to be 
composed of the fuselage unit and the wing unit as in the development 
of the Inertia parameters. On a unit there are acting a normal force n dy 
due to the variation in coordinate Z, a moment m dy about the airplane 
center of gravity due to the variation in coordinate 0 , and a force F dy 
due to the variation in coordinate h. The h coordinate is conposed 
of bending and twisting of the exposed wing. Therefore, the force F dy 
is conposed of force n dy and moment m'dy about the_ local elastic 
axis. The total section work done by the forces on these tuoits corre- 
sponding to a virtual displacement AQ is 


f Q AQ dy = AQ 


3z So 

,n — + m — + 

L Sq Sq 


Sz 

'Sh 


+ m 



(A20) 
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■where 


n = I + 9 ^ + e ^ + h ^ + h ^ 

^ S(z/v) ^ ^0 ^h Sh 


= Z _ ^ +Q ^+0 ^ + 

V S(z/v) Se Se 


m' = 5 _,^J +9^ + e^+’h^ + i^ 
V d(z/v) 50 S 0 ^h 


By definition 


a = - + 0 

V 


Substituting the partial derivati-^s of equation (A2k) into equa- 
tions (A21), (A22), and (A23) yields 


Ti - „ ^ j. a 

n = a + 0 —T + h + h — 
^50 5h 


Sm.T!,.diii 

m = a — + 0 + h — - + h ■:— 

5a 50 5h 5h 


m'=a^ + 0^ + h^ + h^ 

5x 50 5h 5fi 


Upon integrating over the ■wing, equation (A20) becomes 


■p. r'^oM /■ 5 z ^ , 50 \ 5 h 

— + 2/ (n — + m — dv^ 

Sq 5q Jq Bh ^/5 q ° 
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where N and M include the normal-force contrihutions of the fuselsige 
and tail. The integral pertaining to the elastic degree of freedom h 
is limited to the exposed wing span hecause no consonant of the elastic 
motion has heen assumed for the fuselage- Partial differentiation of 
equation (A28) with respect to Z, 6, and h yields 

F2 = N = + 9Nq + hUh + hW^ (A29) 

F0 = M = + 81% + h% + hl^ (A30) 


Fh = 




(A31) 


The generalized force term 
fore^ the coD5)onents of the 

Dividing equation (A31) hy 


Pji is probably the least fam i liar; there- 
force term are obtained in coefficient form. 


q.S and multiplying by 


^o/g 

^0/2 


yields 


qS ICavl'^O N oh/-bj 2 


(A52) 


Let Cp = Fjj/clS, n 
tlon (A32) becomes 


C-nlCj 


t 1 — 


and 


= 


yo 


to/g‘ 


Then equa- 



c 8z 
- 8h 


+ c 


'av 


m 


'av 


_ 80" 




and equations (A25) and (A26) become 


(A35) 


Cq = acn^ + q I Cq^ + Hcn^ + HCn^ (A34) 

Cm' = + q. I + HCm^' + Hcju^* (A 35 ) 
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where by convention, = | Cjn^ ’ = | Cm^' , and H = -^ 

tutlng equations (A5i^) and (A35) into (A35) yields 


Substi- 


Cj. - 




0 L^av 




CC f \ It t • f 1 

^ 15 % 


k„ (»6) 


Partial differentiation of equation (A 36 ) with respect to the 
VEirl OTIS ■ generalized force congjonents a, q_, H, anfl DH, with H = — DH,- 

yields ^ 


Cp 


a 


= — r (cr,-^^+c'-^c I^V^o 
b Jq ^a.^r 3h hhj ° 


(A57) 


I % ' TJ^ (5 % i * 5 %' 


CFh 


■vl (■ 


dZ 


p . I c — 
Crn„ + Cm_ C 


Sh “H Cav bhj 


d’rL (A39) 


%H ' ^ (^^DE ^ S ^ ‘=^’ ° 


(Al^O) 


The equations of motion (eq^s. (AI7) to (AI9) and (A29) to (A31)) 


are now put in nondimensionaJ. form. 


Time is ejqDressed as 


S = 



} 
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the distance traveled in chord lengths. Nondimens ional derivatives are 
obtained as 



Linear quantities are expressed in chord lengths c, and frequency is 
expressed as ta = k ^ ■ Force equations are divided by and moment 

equations by — V^SF. The resulting three equations of motion in non- 
dimensional form are 

2nD(a - 0 ) + - “(%) ' %x) " “(| %) ' 

=(%) - ™(%h) ' 

2|iKy^% - a(b^) - Do(| - I»(| - 

=(%) - ™(Sh) ” 

2AiiD% + 2AziP(a - 0 ) + 2%iiD% + 2*UjA - a(Cj^) - I»(| Cp J - 

h(CFh) - = Cj.^8 (A43) 

The elastic properties of the flexible ••wing mode are manifested 
primarily in the k^ term of equation (A^^5)• The elastic mode shape 
of the flexible -wing is employed in^jlicitly in the generalized mass and 
force terms pertaining to the flexible mode. " 



MCA TN 5251 


25 


'RhiH'kUhilMnK.R 


1. Duncan, ¥. J.: The Representation of Aircraft Wings, Tails and 

Fuselages hy Semi-Rigid Structiires in Dynamic and Static Problems. 
R. & M. No. 1904 , British A.R.C., Feh. 

2. Timoshenko, S.: Vibration Problems in Engineering. Second ed., 

D. Van Rostrand Co., Inc., 1937 • 

3 . Go land, Martin, Luke, Yudell L., and Sacks, Irving: Effects of Air- 

plane Elasticity and Unsteady Flow on Longitudinal Stability. 

Proj. No. R108E 108 (Contract No. AF 33(038)-297^)^ Midwest Res. 
Inst. (Kanseis City, Mo.), Oct. 23, 1950* 

1^. Houholt, John C., and Anderson, Roger A.^ Calculation of Uncoupled 
Modes and Frequencies in Bending or Torsion of Nonuniform Beams. 
MCA TN 1522 , 1948. 

5 . Diederich, Franklin W., and Budiansky, Bernard: Divergence of Swept 

Wings. NACA TN I 68 O, 1948. 

6 . DeYoung, John, and Harper, Charles W. : Theoretical Symmetric Span 

Loading at Subsonic Speeds for Wings Having Arbitrary Plan Form. 
NACA Rep. 921, 1948. 

7 . McLau^lln, Milton D.; Method of Estimating the Stick-Fixed longi- 

tudinal Stability of Wing-Fuselage Configurations Having Unswept or 
Swept Wings. NACA RM L5U23, 1952. 


MCA TN 5251 


TABLE I.- COKFIGDRAnOWS AND PERTINENT CHARACTERISTICS 

[s = l,ll-60 sq^ ft; Z-t = 46.5 ft; e.a. location, 58 percent 
wing chord; wing section c.g. location, 38 percent 
wing chord (approx. J] 


Configuration 


Sweep 

“w/“A 

c.g. 

c, ft 

angle. 

location. 


A, deg 

percent c 


0 

0.15 

25 

11 



35 

11 



45 

11 


.53 

25 

11 



35 

11 



45 

11 


.50 

25 

11 



35 

11 



45 

11 

40 

.33 

25 

14.58 



35 

14.58 



45 

14.58 

60 

.33 

25 

22 



35 

22 



45 

22 i 


Wing ground 
natural frequency, 
CD, radians/sec 


5.98 

5.98 

5.98 


4.21 

4.21 

4.21 


5.98 

5.98 

5.98 


5.98 

5.98 

5.98 





































TABLE n.- nCBZlA A9D DAMPOn AID HKBICEU30-SFHIK] FARAIBEIS 


SVBOP 
angle, 
A, dag 



Farajaatars 


>Mi ^ Sx i %3a I 5 I ^^>*131 


129.5 2.5^ V.6 0.6 3.52 -3.66 -1.23 -1.07 -l^.52 

129.5 2-56 ll-.e .13 3.32 .3.66 -.68 -1,0^ -J1-.3 

129.5 2.36 ^.6 -.32 3.52 -5.66 -.11 -1.02 -4.1 

45 12.53 1.63 2.64 -5.66 -1.25 -1.07 -4.52 

45 12.55 .3S 2.64 -3.66 -.66 -1.04 -4.3 

43 12.33 -.06 2.64 -3.66 -.U -1.02 -4.1 

9 23.2 3-27 1.96 -5.66 -1.25 -1.07 -4.52 

9 25.2 .75 1.98 -3.66 -.68 -l.o4 -4.3 

9 25.2 -1.76 1.98 -3.66 -.11 -1,02 -4.1 

125.0 4.93 9.61 9.61 1.75 -4.24 -.93 -.82 -2.63 

125.8 4.93 9.61 8.63 1.73 -4.24 -.51 -.79 -2.49 

123.0 4.93 9.61 7.68 1.75 -4.£^^ -.08 -.77 -2.33 

22 6.28 6.07 -2.98 -.65 -.55 -1.13 

22 6.28 3.44 .85 -2.98 -.56 -.31 -1.02 

22 6,20 4.82 .63 -2.90 -.06 -.48 -.92 


-3.96 -11.69 0.042 0 -1.39 0 
-3.44 -11.17 .042 0 -1.39 0 
-2.93 -10.65 .042 0 -1.39 0 

-3.96 -11.69 .042 0 -1.39 0 
-3.44 -11.17 .042 0 -1.39 0 
-2.95 -10.63 .042 0 -1.39 0 

-3.96 -11.69 .042 0 -1.39 0 
-3.44 -11.17 .042 0 -1.39 0 
-2.93 -10.63 .042 0 -1,39 0 


-.40 -.268 

-1.00 

-.95 

-.40 -.268 

-1.00 

-.95 

-.40 -.268 

-1.00 

-.95 

-.30 -.519 

-.62 

-.564 

-.50 -.519 

-.62 

-564 

-.50 -.319 

-.62 

-364 


-1-39 0 
-1.39 0 
-1.39 0 

-1-39 0 
-1-39 0 
-1.39 0 

-1.59 0 
-1.39 0 
-1.39 0 

- 1.00 0 
- 1.00 0 
- 1.00 0 

-.62 0 
-.62 0 
-.62 0 


i Cp 



-0.224 

0.0:9 

“0.77 

-.134 

.019 

-77 

-.043 

.019 

-77 

-.224 

•ois 

-.77 

-.134 

.019 

-.77 

-.045 

.019 

-77 

-.224 

.019 

-•n 

-.134 

.019 

-.77 

-.045 

.019 

-77 

-.848 

n.177 

-35 

-.705 

-.177 

-.35 

-.719 

-.177 


-.416 

-.222 

-.34 

-.578 


-.34 

-.342 

-.222 

-34 
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Figure 1 .- A diagram showing the stahlUty axes and positive directions 

of forces and moments. 
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0 ,2 ,ij. .6 



0 .2 .U 


Spanwlse station, 

Figure 3 -- Wing spanwise elastic and weight distributions 
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(a) Frequency and damping of elastic wing mode and frequency of airplane 
short-period longitudinal stahility mode in dimensional terms as a 
function of dynamic pressure. 

Figure 5 .- Eynamlc stability of airplane longitudinal mode and flexible 
wing mode obtained by the semirigid method and less rigorous methods 
as a function of dynamic pressure. Altitude, 8,000 feet; A « 0° . 








I%rlod^ Inswltge langthfl 



test VDVH 










Tlmd to 0.1 amplitudOf Tq boc Freituoacyi ops 


34 


MCA TN 3251 



0 200 <400 600 800 1000 0 200 400 600 800 1000 


Dynamic pressare, q, U/aqfi Dynamic p j e aaur s, q»lb/sqlt 



0 200 400 000 800 1000 0 200 400 600 800 1000 


Dynamic pressure, Dy nandc pr e Bsnr e, q,lV&4ft 

(a) Frequency and dar^jing of elaistic wing mode in dimensional terms 
ELS a function of dynamic pressure. 

Figure 6 .- Dynamic stability of airplane longitudinal mode and flexible 
id.ng mode obtained by the semirigid method and less rigorous methods 
as a function of dynamic pressiore. Altitude, 8,000 feet: A = 4o° 
and 60 °. 
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Maneuver margin, chords 


deg 

m„/m^ 

c.g. 

position 

%c 

0 

.33 

25 

0 

.33 

35 

60 

.33 

25 

ho 

.33 

25 

0 

.33 

li5 

60 

.33 

35 

ho 

.33 

35 

60 

.33 

h5 

ho 

.33 

h5 


,00 200 300 UOO 500 600 700 800 900 

Dynamic pressure, q, Ib/sq ft 

Flgirre 7.- Variation of the maneuver margin with dynamic pressure for 
configurations with flexible wings. Quasi-static method. Altitude, 
8,000 feet. 
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(a) Altitude, 8,000 feet; Mach number, 0.4; q = I 76 Ib/sq ft. 

Figure 8.- Longitudinal transient responses of typical configuration 
with flexible wing to 1-radian step inputs of the elevator at 
altitudes of 8,000 feet and 30,000 feet. Center-of -gravity location, 
20 percent mean aerodynamic chord; A = 55° • 



Angle of Normal 

Wing- tip deflection, attack, acceleration, Pitching velocity, 

h, ft a, radians g units 0, radians/sec 
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0 123J+5678 

Time, t, sec 


(b) Altitude, 30,000 feetj Mach number, 0.6j q. = 158 Ib/sq ft. 

Figure 8.- Concluded. 
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(a) GeometirLc relation of fuselage center of gravity to airplane 

center of gravity. 



z„ f^(y)h z 


(b) Geometric relation of local wing center of gravity to airplane 

center of gravity. 

Figure 9 *- Linear and angular velocities and geometric relations of the 
fuselage and wing centers of gravity to the airplane center of gravity. 
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